Microscopic Dynamics in a Strongly Interacting Bose-Einstein Condensate 
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An initially stable 85 Rb Bose-Einstein condensate (BEC) was subjected to a carefully controlled 
magnetic field pulse in the vicinity of a Feshbach resonance. This pulse probed the strongly inter- 
acting regime for the condensate, with calculated values for the diluteness parameter (na 3 ) ranging 
from 0.01 to 0.5. The field pulse was observed to cause loss of atoms from the condensate on re- 
markably short time scales (>10 fj,s). The dependence of this loss on magnetic field pulse shape and 
amplitude was measured. For triangular pulses shorter than 1 ms, decreasing the pulse length actu- 
ally increased the loss, until extremely short time scales (a few tens of microseconds) were reached. 
Such time scales and dependencies are very different from those expected in traditional condensate 
inelastic loss processes, suggesting the presence of new microscopic BEC physics. 

PACS numbers: 03.75.Fi, 05.30.Jp, 32.80.Pj, 34.50.-s 



The self-interaction strength that determines most of 
the properties of a Bose-Einstein condensate (BEC) is 
characterized simply by the s-wave scattering length, a 
H . Near a Feshbach resonance Q , the scattering length 
depends strongly on the magnetic field. This property 
has been used to vary the self-interaction energy in BEC 
H |J, H ||, 0- In particular, we have used the Feshbach 
resonance to create stable 85 Rb condensates ||. We also 
varied the scattering length and studied the effect on the 
condensate — most notably its collapse when the scat- 
tering length was made negative |j, [tJ . 

Here we discuss BEC behavior in the positive a re- 
gion near the 85 Rb Feshbach resonance, where the self- 
interaction is large and repulsive. We used rapid mag- 
netic field variations to probe the strongly interacting 
regime in the condensate and to investigate the possi- 
bility of collisional coupling between pairs of free atoms 
and bound molecules [§, §, [l| H H- We changed 
the magnetic field to approach the Feshbach resonance 
from above, which causes the positive scattering length 
to increase. As na 3 becomes comparable to or larger 
than one, the customary mean-field description of the di- 
lute gas BEC, which assumes no correlations between the 
atoms, breaks down. Correlations between the atoms be- 
come increasingly important and a BEC with sufficiently 
large na 3 will eventually reach a highly correlated state 
such as that found in a liquid. We are interested in the 
nature of this transition to a highly correlated state and 
the time scale for the formation of the condensate cor- 
relations. In addition to such interesting physics, our 
time dependent experiments also allow us to probe the 
coupling between atomic and molecular BEC states that 
are degenerate at the Feshbach resonance. Several au- 
thors have proposed mechanisms for transferring atom 
pairs into molecules using time-dependent magnetic fields 
near resonance [|[ [| |l(| [n], [l^, [lj| . To investigate such 
physics we examined the response of the condensate when 
we briefly approached the resonance by applying a short 



magnetic field pulse (much shorter than the period for 
collective excitations) . The most obvious feature of the 
BEC response was a loss of atoms that increased when 
we made our pulses shorter, until the loss became very 
large for very short pulses. 

Bose-Einstein condensates always disappear if one 
waits long enough. The dominant loss process was found 
to be three-body recombination into molecules Jl4], [TBI , 
with a time dependence characterized by a simple rate 
equation and rate constant. All previous observations 
of condensates are consistent with a mean-field descrip- 
tion that includes such a density-dependent loss process. 
In the vicinity of a Feshbach resonance the decay rates 
have been seen to increase dramatically ||, |[ |lq| , be- 
coming so large that various novel coherent conversion 
processes have been put forth to explain their magnitude 
@, 0, |§. However, the observations in Refs. |, |, |l|] 
always revealed that the more time spent near the Fes- 
hbach resonance, the greater the loss. This is still con- 
sistent with the basic picture of a mean-field with some 
inelastic loss term, albeit a very large one. Here we re- 
port losses occurring on very fast time scales even when 
we remain some distance from the resonance, and we see 
that shorter and more rapid pulses lead to a greater loss 
than longer, slower pulses. This behavior contrasts dra- 
matically with the conventional picture presented above, 
and indicates the presence of novel BEC physics. 

To study BEC dynamics we first formed a 85 Rb con- 
densate, following the procedure given in Ref. ||. A sam- 
ple of 85 Rb atoms in the F=2 mp=-2 state was evapora- 
tively cooled in a cylindrically symmetric magnetic trap 



dial 



17.5 Hz 



6.8 Hz). The magnetic 



field at completion of evaporation was 162.3 G, corre- 
sponding to a scattering length of 200 ao- Typically, the 
cooled sample had A*o=16,500 BEC atoms and fewer than 
1000 noncondensed atoms. The magnetic field was then 
ramped adiabatically (in 800 ms) to ~166 G where the 
scattering length was positive but nearly zero |17|, and 
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the BEC assumed the shape of the harmonic trap ground 
state. 

We next applied a short magnetic field pulse (dura- 
tion < 1 ms) so that the field briefly approached a value 
moderately close to the Feshbach resonance at ~155 G 
(see Fig. 1). We used destructive absorption imaging to 
look at the number of atoms remaining in the conden- 
sate. This experiment was repeated with a variety of 
differently shaped magnetic field pulses. 

We found that the magnetic trap must be turned off 
and the condensate spatial size must be significantly 
larger than our resolution limit to obtain the most sen- 
sitive and accurate measurements of number [Q . To ex- 
pand the BEC after the short pulse, we ramped in 5 ms 
to ~157 G (a=1900 a ), and then held at that field for 
7 ms. The mean-field repulsion during the ramp and hold 
times decreased the density by about a factor of 30, then 
we rapidly turned off the trap and imaged the cloud. The 
density decrease avoided density-dependent loss that we 
have observed during the trap turn-off j|. [Hi] . 

Significant number loss from the BEC was observed for 
pulses lasting only a few tens of microseconds. The loss 
was accompanied by the generation of a "burst" of rela- 
tively hot (~150 nK) atoms that remained in the trap, 
as in Ref . . The few thousand burst atoms represented 
a fairly small fraction of the total number lost from the 
BEC and so in this Letter we have focused only on the 
number in the condensate remnant. The burst will be 
the subject of future work. 

To study the remarkably short time scales for the loss, 
we designed a low-inductance, high current auxiliary elec- 
tromagnetic coil. The coil current was supplied by a ca- 
pacitor bank that was charged to 580 V, then discharged 
through the coil at a rate controlled by a transistor. Our 
goal was to create a perfect trapezoidal magnetic field 
pulse with adjustable but identical rise and fall times 
and a hold time during which the field was constant. 
To compensate for observed mutual inductance effects 
with other coils in the apparatus, we empirically deter- 
mined the auxiliary current pulse shape needed to pro- 
duce a total magnetic field pulse that closely approxi- 
mated our ideal, as shown in Fig. 1. The presence of 
induced currents limited the maximum ramp speed to 
dB/dt =1 G//is, and the magnetic field uncertainty was 
1 part in 10 3 (0.16 G). 

We first examined the BEC loss for trapezoidal field 
pulses as a function of the hold time (see Fig. 2). Using 
a linear rise and fall time of 12.5 /is, we observed BEC 
number loss of 10-20% when the hold time was set to 
zero (triangular pulse). The number of BEC atoms then 
showed a smooth exponential decrease as hold time was 
increased. Surprisingly, when we reduced the initial BEC 
density by more than a factor of 2, the time constant was 
nearly unchanged. For the low[high] density data, the 
value of the time constant was over 2[1] orders of mag- 
nitude shorter than predicted by our previous inelastic 
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FIG. 1: The magnetic field vs. time for a typical auxiliary 
coil pulse. The solid line is the field calculated from mea- 
surements with calibrated Hall-effect current sensors on each 
electromagnet coil. On the right vertical axis, the correspond- 
ing variation in scattering length is shown. The open circles 
are independent checks of the magnetic field obtained from 
the resonant frequency for a 10 /is RF radiation pulse that 
drives atoms to the niF=-l spin state (error bars are smaller 
than the points). The dashed line shows the position of the 
Feshbach resonance, where the scattering length becomes infi- 
nite. The field variation on the peak of the pulse was typically 
AB < 0.1 G. 



decay measurements with cold thermal clouds |18 



We next measured how the loss depended on the 
rise/fall time of the pulse for a variety of pulse ampli- 
tudes and hold times (Figs. 3 and 4). We varied the 
rise time from 12.5 /is to 250 fis and changed the hold 
time at the pulse peak from 1 fis to 100 /is. In addition, 
the amplitude was varied to examine fields from 158.0 G 
(a=1100 a ) to 156.0 G (a=4000 a ). Here we list the 
corresponding scattering lengths that we have observed 
by slowly adjusting the magnetic field, as in Ref. [[|. 
The scattering length was calculated from the equation: 
a(B) = a bg (l - A/{B - B )), where A=11.0(4) G is 
the width of the resonance, i?o=T54.9(4) G is the res- 
onant magnetic field, and a;, g =-450(3) ao is the back- 
ground scattering length |Q. For the range of magnetic 
fields examined here, the initial value of the dilutcness 
parameter varied from na 3 = 0.01 for a=1100 ao to na 3 
= 0.5 for a=4000 a . 

The number remaining in the BEC after the pulse 
vs. pulse rise time for a variety of different hold times is 
shown in Fig. 3. For hold times thoid < 15 /US, there is an 
initial decrease in N rem as the rise time increases. Then 
the slope changes and fewer atoms are lost for longer 
rise times. All of the hold time data display this up- 
ward slope over some range, but the range is largest for 
the 100 /is hold time. The increase in remnant number 
for longer rise time provides clear evidence that the loss 
is not conventional inelastic decay that is characterized 
by a rate constant. It is interesting that the short hold 
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FIG. 2: Fraction of BEC remaining vs. pulse hold time. The 
pulse rise/fall time was 12.5 fis and the magnetic field during 
the hold was 156.7 G (2300 ao). Number decay was mea- 
sured for two different initial densities: (n) = 1.9xl0 13 cm -3 
(•) and (n)=0.7xl0 13 cm -3 (o). Fitting the data to exponen- 
tial functions (solid lines) gave time constants of 13.2(4) /is 
and 15.4(14) /is. Thus, reducing the density by a factor of 2.6 
caused an increase of only 17(11)% in the time constant. 



time data show a distinct minimum in N rem vs. rise time, 
which shifts toward shorter rise times as the hold time is 
increased. 

In Fig. 4 we display the remnant number vs. ramp 
time for triangular magnetic field pulses (1 /is hold). The 
remnant number is shown for various pulse amplitudes. 
For all cases, N rem is large at the shortest rise time and 
then decreases with rise time until it reaches a minimum. 
Then longer rise times cause N rem to increase over a time 
scale of tens of microseconds. The rise time that induces 
maximal loss becomes longer as the pulses come closer to 
resonance. 

Conventional condensate loss is characterized by a rate 
constant for a density-dependent decay process, and thus 
the loss increases monotonically with time. Near a Fesh- 
bach resonance the rate constants have been observed to 
increase enormously ||, ||, but nevertheless, a longer 
time spent near the peak, or equivalcntly a slower ramp 
getting over it, resulted in greater loss. In contrast, we 
have measured an increase in the loss when the ramp 
time is decreased, which reveals the existence of previ- 
ously unexplored BEC physics. Of course, the above in- 
terpretation of our data would be modified if the density 
of the sample was changing due to the rapid increase 
in the mean-field interaction, but the characteristic time 
for such readjustments in cloud shape is far longer (of 
order 1/ ' {2v ra di a {)=2 < d ms) than the time scales we con- 
sidered here. For example, using a simple analytic model 
based on the Gross-Pitaevskii equation p3], we calcu- 
late that for a 250 /is ramp to Bfinai^ 156.0 G (4000 ao), 
the change in mean-field energy causes the density to 
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FIG. 3: Dependence of remnant BEC number on pulse 
rise/fall time for various hold times (see legend) with 
iV =16,500 ((n) =1.9xl0 13 cm" 3 ). For the majority of the 
data points, the symbol is larger than the statistical error 
bar (not shown). The lines are spline fits to guide the eye. 
The abscissa was multiplied by a factor of 1/4 to show the 
time required to ramp from 75% to 100% of the pulse ampli- 
tude. This scaling reflects the observed fact that most of the 
loss occurred at fields closest to the Feshbach resonance. The 
magnetic field during the hold time was 156.7 G (2300 ao). 
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FIG. 4: Rise time dependence of BEC number for different 
triangular pulse amplitudes. The pulse hold time was fixed 
at 1 /is for these data and No = 16,600. Data points are larger 
than the error bars. 



decrease by only 1% from its initial value. Thus, the ob- 
served dependencies on rise time must reflect microscopic 
physics in the BEC and not any macroscopic changes in 
the shape of the condensate. 

The response of the BEC remnant to magnetic field 
ramps, observed in Fig. 4, is qualitatively similar to what 
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one would expect in a Landau-Zener (L-Z) model |g3|, |24[] 
of an avoided crossing of two states with linear Zeeman 
shifts, if the second state were undetectable. In our case, 
the obvious candidate for this second state is the Fesh- 
bach resonance bound state corresponding to a diatomic 
molecule. With a L-Z avoided crossing, the behavior 
when the field approaches and then backs away from the 
crossing point with a triangular pulse is qualitatively sim- 
ilar to the more familiar and analytically soluble case ]2j| 
of a linear field ramp that goes from far above the cross- 
ing to far below and then back again. In both cases, 
and as seen in Fig. 4, the L-Z model predicts a steep rise 
from zero in the transition probability as the length of the 
ramp increases from zero (diabatic limit). As the ramp 
time increases further, the transition probability reaches 
a maximum — where the time derivative of the relative 
energy matches the square of the coupling strength - 
and then slowly decreases to zero as the ramp approaches 
the adiabatic limit. We found that the L-Z model does 
a rather bad job at reproducing any more quantitative 
features of the data however, even when we allowed the 
coupling strength and relative magnetic moment to be 
arbitrary parameters in our numerical simulations p(| . 

Nevertheless, it seems likely that some atoms are being 
converted to another state (possibly molecular) by nona- 
diabatic mixing, although the process is more compli- 
cated than a simple Landau-Zener avoided crossing pic- 
ture. In future experiments, we plan to further investi- 
gate the time response of the BEC loss using asymmetric 
pulses and double pulses with variable spacing. We will 
explore the burst production process and attempt to de- 
termine the fate of the lost atoms. 
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